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The early dynamical evolution of cool, clumpy star clusters 
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ABSTRACT 

Observations and theory both suggest that star clusters form sub-virial (cool) with 
highly sub-structured distributions. We perforin a large ensemble of iV-body simu- 
lations of moderate-sized {N — 1000) cool, fractal clusters to investigate their early 
dynamical evolution. We find that cool, clumpy clusters dynamically mass segregate 
on a short timescale, that Trapezium-like massive higher-order multiples are com- 
monly formed, and that massive stars are often ejected from clusters with velocities 
> 10 km s~^ (c.f. the average escape velocity of 2.5 km s~^). The properties of clus- 
ters also change rapidly on very short timescales. Young clusters may also undergo 
core collapse events, in which a dense core containing massive stars is hardened due 
to energy losses to a halo of lower-mass stars. Such events can blow young clusters 
apart with no need for gas expulsion. The warmer and less substructured a cluster is 
initially, the less extreme its evolution. 
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1 INTRODUCTION 



Most sta rs appear to for m in star c lusters ("Lada fc Ladal 
l2003l : [Lad a 2009; Portegi es Zwart et a l. 2010) and so star 
formation is inextricably linked with star cluster formation. 
Recent advances in observations and theory have allowed us 
to construct a basic picture of cluster formation in which 
clusters form dynamically cool (sub-virial), and highly sub- 
structured. 

Clusters form in highly turbulent molecular clouds. 
These clouds are hig hly substructured, containing dense 
clumps and filaments (IWilUamsHigggl : IWiUiams et aLlbOOd: 
ICarpenter fc Hodapd [2008) which are presumably formed 
by the decay of supersonic turbulenc e (iMac Low fc KlessenI 
l2004l : iBallesteros-Paredes et al.ll2007h . Stars and small stel- 
lar groups form in these dense regions and are unsur- 
prisingly ob served to haye a high degree of substructure 



lOffner et al] 120091 ') . including in comp arison tests between 
AMR and SPH techniques (|Federrath et al..,20ld '). 



Clusters are observed to lose their substructure as 
they evolve, becoming s moot h and roughly spherica l 
feartwright fc WhitworthI |2004|: ISchmeia fc KlessenI l2006l : 



iSchmeia et al.ll2008l ). ICoodwin fc WhitworthI (|2004r ) have 
shown that substructure can only be erased in clus- 
ters if the clusters are initially cool (sub-virial) (see 
also Maschb erger et al. 2010). Both observations of pre- 
stellar cores (iBelloche et al.ll200ll: lAndrell2002l: fWalsh et al. 



when young jLarso n 1995*: Elmcgrccn '200^; Testi et al 



Cartwright fc Whitworth ,2004: ,Gutermuth et al 



200C 



2005 



Allen et al.l l2007l : ISchmeia fc KlessenI l2006l : ISchmeia et al.l 
20081 '). clustering in you ng stellar groups ha s even been 



observed in the SM C l|Schmeia et al.l 120091 ') and LMC 
jBastian et al.ll2009l ). The same behaviour is seen in sim- 
ulations of cluster formation (iKlessen fc BurkertI 2000l . 
I2OOII : iBate et al.l l2003l : iBonnell eTal.! l2003l : iBatd l2009l : 
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I2OO4I: jPeretto et al.l l2006l: iKirk et al.l |2007^ and stars 
jPeretto et al.ll2006l : [Proszkow et ablbOOgf show that they 
indeed appear to be sub-virial, a property also found in 
simulations of cluster formation (Klessen fc Burkert 2000; 
Offner et al. 2009; Maschberger et al. 2010). 



Following Allison et al. (2009b) we conduct TV-body 
simulations of a large number of initially sub-virial, frac- 
tal star clusters. In Section [2] we describe our simulations. 
In Section [3] we describe our main results, specifically the 
early onset of dynamical mass segregation and interesting 
'post-collapse' evolution. In Section|4]we discuss the results, 
and we summarise and conclude in Section [5| 
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2 INITIAL CONDITIONS 

We perform 160 A''-body simulations with 1000 stars each, 
in which the initial conditions are cool and clumpy. We vary 
the level of substructure and initial virial ratio, and conduct 
ensembles of simulations with the same initial conditions, 
varying only the initial random immber seed used to ini- 
tialise the simulations. 

To create initial substructure in our simulations we 
use a fractal stellar distribution. Using a fractal distri- 
bution provides a parameterisation of substructure using 
only a single number: the fractal dimension. (Note that we 
are not claiming that clust ers are actually initially frac- 
tal, although they may be (lElmegreen fc ElmegreenI l200ll : 



ICartwright fc WhitworthI 12004 ). just that this provides a 
simple descriptor of substructure that is easy to reproduce). 

The fractal stella r distributions were generate d fol- 
lowing the method of iGoodwin fc WhitworthI (|2004l ). The 
method begins by defining a cube of side A'div (we use 
A^div = 2 throughout), inside of which the fractal will be 
built. A first-generation parent is placed at the centre of the 
cube, from which are spawned N^i^ sub-cubes, each con- 
taining a first-generation child in its centre. The fractal is 
then built by determining which of the children themselves 
become parents, and spawn their own offspring. This is de- 
termined by the fractal dimension, D, where the probabil- 
ity that a child becomes a parent is Nj^^~^K For a lower 
fractal dimension less children will mature and so the final 
distribution will contain more structure. Any children which 
do not become parents in a given step are removed, along 
with all of their parents. A small amount of noise is then 
added to the positions of the remaining children, prevent- 
ing the final cluster from having a gridded appearance, and 
the children become parents of the next generation. Each 
new parent then spawns A^l;^ second-generation children in 
A'div sub-sub-cubes, with each second-generation child hav- 
ing a N^^~^^ probability of becoming a second-generation 
parent. This process is then repeated until there are substan- 
tially more children than required. The children are pruned 
to produce a sphere from the cube and are then randomly 
removed (so maintaining the fractal dimension) until the 
required number of children are left. These children then 
become the stars in the cluster. 

To determine the velocity structure of the cloud, chil- 
dren inherit their parent's velocity plus a random component 
that decreases with each generation in the fractal. The chil- 
dren of the first generation are given random velocities from 
a Gaussian of mean zerc0. Each new generation then inher- 
its their parent's velocity plus an extra random component 
that becomes smaller with each generation. This results in 
a velocity structure in which nearby stars have similar ve- 
locities, but distant stars can have very different velocities. 
Finally, the velocity of every star is scaled to obtain the 
desired total virial ratio for the cluster. 



The simulations contain 1000 stars, have an initial max- 



^ Here the variance is unimportant as the velocities are scaled 
to the desired virial ratio once the final spatial distrib ution is 
obtained, but the method can also be used to match the [Larson] 
lligSlI) relations if desired fsee lGoodwin fc Whitworthll2004l ') . 
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Table 1. Notation for run identification where D is the initial 
fractal dimension, and Q is the initial virial ratio of each simu- 
lation. The numbers 01—50 and 01-10 are the identifiers for each 
individual run. Within each ensemble only the random number 
seed used to generate the initial conditions is changed. 



imum radius of 1 pc, include no primordial binaries or gas 
and a three-par t power law is used to produce an initial mass 
function (IMF. lKroupall2002l '). 



r M-°-^ mo ^ M/ M0 < mi, 
N{M) oc <^ mi ^ M/ M0 < ma, 



(1) 



1712 ^ M/ Mq < Tn3, 



with mo = 0.08 Mq, mi = 0.1 Mq, ma = 0.5 Mq and 
ma = 50 Mq. No stellar evolution is included because of 
the short duration of the simulations (~ 4 Myr). We use 
the STARL AB A^-body integrator kira to run our simulations 
(Portcgic s Zwart et al, 2001) . 

In this study we explore a range of fractal dimensions 
and virial ratios. The fractal dimensions investigated are 
D — 1.6, 2.0, 2.6 and 3.0 (since these values correspond 
to the number of maturing children, A^^iv = 2^ , being an 
integer), where D — 1.6 produces a large amount of struc- 
ture, and D = 3.0 produces a uniform sphere. We investigate 
virial ratios of Q = 0.3, 0.4 and 0.5, we define the virial ra- 
tio as Q = T/\^l\ (where T and \Q\ are the total kinetic and 
total potential energy of the stars, respectively), hence virial 
equilibrium is Q = 0.5. 

It is important to note that fractal initial conditions are 
inherently stochastic: statistically identical fractals (i.e., the 
same fractal dimension), can appear very different to the 
eye, and can evolve in very different ways (see Section [3. 3p . 
Therefore, it is vital to perform large ensembles of simula- 
tions with different random number seeds. We have therefore 
simulated 50 D = 1.6, Q = 0.3 (al) clusters (as they have 
the most interesting evolution), and restricted our analysis 
of all other combinations of D and Q to 10 clusters each. 
The identifiers and initial conditions of each ensemble are 
presented in table [T] 

To quantify mass segreg atio n we use the minimum span- 
ning tree (MST) method C Allison e t al, 2009a). This method 
compares the MST of the A'^ most massive stars with the av- 
erage MSTs of A'^ randomly selected stars. The ratio of these 
two MST lengths gives a quantitative measure of the con- 
centration of the massive stars and hence a value of the mass 
segregation in the cluster, 



A 



(^norm) 



massive 



(2) 



massive 



where A is the measure of mass segregation, (Znorm) is the 
average length of the random MSTs, /massive is the length of 
the massive star MST and (Tnorm is the error in the random 
MST length. The value cTnorm/^massive is the 1 a error in A. 

In this paper, mass segregation is calculated for the en- 
tire stellar content (i.e. ejected stars are not removed from 
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the simulation) in all three dimensions. We note that an 
observer will see only two dimensions, and may not iden- 
tify all cluster members, in particular low-mass stars and 
stars that have been ejected from the cluster. Allison et al. 
(2009a) show that there is little difference between the mass 
segregation measure in two and three dimensions when us- 
ing mass segregated Plummer spheres. However, as with all 
mass segregation measures, if using the measure in two di- 
mensions projection effects may cause differences depending 
on orientation. Unfortunately, this is unavoidable. Including 
all stars is useful theoretically, however it does not allow a 
direct comparison with observational data. 



3 RESULTS 

3.1 Rapid dynamical mass segregation 

InlM ison et al. I l|2009bl . hereafter Paper I) we showed that 
cool (Q = 0.3), and clumpy [D = 1.6) clusters can dynami- 
cally mass segregate on a timescale close to the initial cross- 
ing time of the system (~ 1 Myr). Figure [T] shows the early 
evolution and onset of mass segregation in such a cluster. 
The top panels show the spatial distributions of the cluster 
at times ~ 0.01,0.81 and 1.02 Myr (left to right), whilst 
the lower panels show the mass segregation ratio A for each 
of these snapshots - the higher the value of A, the greater 
the degree of mass segregation. As is clear from the figure, 
the cluster evolves from an initially non-mass segregated and 
clumpy distribution, into a smooth and mass segregated one. 

In paper I we argued that this rapid dynamical mass 
segregation is due to the collapse of the cool cluster form- 
ing a very dense, but short-lived core (at ~ 0.8 Myr in the 
example shown in Fig. [T|. These dense cores (for D — 1.6, 
Q = 0.3) tend to contain about half of the mass of the clus- 
ter, are roughly 0.1 - 0.2 pc in radius, and survive for 0.1 - 
0.2 Myr (10 - 20 c rossing times in the core). 

ISpitzed (Il969l ) showed that the timescale for mass seg- 
regation, iaeg, for a star of mass M depends on how massive 
that star is relative to the average mass of a star in a cluster, 
(m). 



tseg{M) 



(m) 
M 



trcla 



(3) 



where trciax is the relaxation time of the cluster, which 
is related to its crossing time, tcross, by 

N 



^ ' 81niV*"°==- 
Eq. 13] can be rewritten as 

7 ''C 



(4) 



(5) 



For the D = 1.6, Q = 0.3 initial conditions used in 
paper I, typical values for these parameters are A'^ ~ 300 
- 500, tcross ^ 0.01 - 0.2 Myr, (m) ~ 0.4 Mq (typical for 
standard IMFs). The core has a lifetime of 0.1 - 0.2 Myr in 
which it can mass segregate giving a mass to which the core 
can segregate of M ~ 2 - 4 Mq. In Fig.[T]the cluster is mass 
segregated to the 20th most massive star, which has a mass 
of a few Mq, in good agreement with our simple analytical 
model. 

The reason that a dense core can form is that the cluster 
is both cool and clumpy: cool clusters will initially collapse, 



but cool and clumpy clusters can collapse further. The po- 
tential energy, il, of a cluster of mass A/cius and radius R 
is 



GM. 



clus 



R 



(6) 



where t; is a structure parameter whose value depends 
on the choice of R (e.g., is it the core radius or the 
half- mass radius?), and the structure of the cluster (e.g., 
is it clumpy, Plummer, o r uniform density?) (see also, 
IPortegies Zwart et al.l[20ld ). For example, for a Plummer 
sphere, if R is the Plummer radius then rj ~ 0.3. 

If a cluster has an initial potential energy flo (with ra- 
dius _Ro and structure parameter -qo), and an initial virial 
ratio Qo, then the initial total energy Eo is 



Eo = -Vo 



Ro 



(1-Qo). 



Whatever the initial conditions, a (bound) cluster will 
attempt to reach virial equilibrium and a relaxed configura- 
tion (something like a Plummer sphere, or King profile with 
a concentration parametei0 ~ 2 — 3, for an al type cluster). 
Therefore, the final energy -B/ of the cluster (with potential 
energy Qf, radius Rf, structure parameter rjf, and virial 
ratio Qf — 0.5) will be 



1 GMl^, 



Rf 



(assuming no mass is lost). Equating these equations gives 
the degree of collapse (or expansion if initially warm) of the 
cluster 



(7) 



Clearly, to induce dynamical mass segregation on a 
short timescale, Ro/Rf must be as large as possible to cause 
the maximum degree of collapse. This implies that both 
rjo/rif and 2(1 — Qo) need to be large. 

Obviously, making 2(1 — Qo) large means making Qo 
as small as possible, but even with an initially static stellar 
distribution with Qo — 0, 2(1 — Qo) can never be greater 
than two (this is why Bonnell & Davis 1998 did not see 
rapid, early mass segregation). Realistically, it is difficult to 
imagine Qo being much less than 0.2 or 0.3 as the stars must 
initially have some relative velocities. 

As stated above, the cluster will not just relax into virial 
equilibrium, it will also attempt to reach a basic statistical 
equilibrium which is a smooth, centrally concentrated dis- 
tribution (like a Plummer sphere or King model if tidally 
truncated). Therefore the final structure parameter will be 
rjf ~ 0.3 when the scale radius is the Plummer radius. 

Numerical experiments we have carried out show that 
the structure parameter for a fractal of dimension D = 1.6 
is ?7o ~ 1.1 ± 0.1, for D = 2.0, 770 ~ 0.8 ± 0.1 (note that 
different realisations of D = 1.6 and D = 2.0 fractals have a 
large variation in their structure parameters), for D = 2.6, 
rjo ~ 0.7, and for D = 3.0, rjo ~ 0.6. For each of these Tj 
values the scale radius is the total radius of the fractal. Un- 
fortunately, it is very difficult to exactly compare the 77 of 

^ We define the concentration parameter as c = 

loglo('"virial/»'corc) 
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Figure 1. Run al.34 - Top: 2 dimensional stellar distributions at (a) Myrs - the initial distribution, (b) 0.8 Myrs - the formation of 
the dense core and (c) 1 Myr - the end of the dense core and initial dynamical mass segregation. Triangles indicate the positions of the 
10 most massive stars. Bottom: The evolution of A with N for (a), (b), (c); as described above. The error bars show Icr deviations. 



a Plummer model and the r; of a fractal as the radii are 
defined differently. The half-mass radius of fractal models 
varies quite significantly (and its meaning is also rather un- 
clear in a fractal) and so does not provide a useful compar- 
ison radius either. 



3.2 Mass segregation in clusters with different D 
and Q 

The analysis above suggests that the initial virial ratio and 
the initial fractal dimension of a cluster are the crucial pa- 
rameters in determining whether that cluster will be able 
to rapidly dynamically mass segregate. In particular, the 
warmer (higher-Q), and the smoother (higher-D) a cluster 
is, the lower the maximum central density, and hence less 
dynamical mass segregation will occur. 

In Fig.[2]we show the typical evolution of A with time for 
clusters with different initial virial ratios and fractal dimen- 
sions for A'' — 10, 20 and 50, the full version of this figure, 
showing the evolution of all of the simulations, can be found 
in the supplementary data. In the top left is our canonical 
Q = 0.3, D — 1.6 cluster from paper I. From left to right, 
the initial virial ratio increases from Q = 0.3 to 0.4 and 
0.5 (virialised) . From top to bottom, the fractal dimension 
increases from D = 1.6 (very clumpy) to 2.0, 2.6 and 3.0 
(roughly a uniform density sphere). 

It would be expected from our earlier argument that 
lower-Q and lower-D clusters will collapse to a denser state, 
and hence show more rapid and more pronounced dynamical 



mass segregation. This is exactly what is seen in Fig. (2] faster 
and more intense mass segregation at the top left (cool and 
clumpy clusters), and no appreciable mass segregation at 
all at the bottom right (virialised, uniform density clusters). 
The trend in the typical evolution of clusters in our param- 
eter space shown in Fig. [2] is exactly what we would expect 
to see following the theoretical argument in Section [3TT] 

Whilst the behaviour of 'typical' clusters is exactly what 
is expected, many individual clusters show unusual and un- 
expected behaviour. Examination of each simulation (see the 
supplementary data) shows that, whilst the general pattern 
of evolution with D and Q holds, there is a large degree of 
stochasticity due to each fractal being different in detail to 
every other fractal. 



3.3 Stochasticity 

The supplementary data shows that a number of simulations 
do not show the behaviour that we might expect for their 
D and Q values. The evolution of fractal clusters depends 
significantly on the specific initial conditions of the cluster. 
The systems are inherently chaotic and so each new random 
number seed can deliver completely different evolution. Phe- 
nomena which occur during the evolution of a cluster (such 
as mass segregation) can also be very transient. In some sys- 
tems, A may change by a large factor on short timescales. In 
this section we have chosen a few examples which demon- 
strate the stochasticity and transience in the evolution of 
cool, clumpy clusters. 
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(a) D = 1.6; Q = 0.3 - al.21 



Time (Myrs) 

(b) D = 1.6; Q = 0.4 - bl.05 
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(d) D = 2.0; Q = 0.3 - a2.08 
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(e) D = 2.0; Q = 0.4 - b2.08 
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(c) D = 1.6; Q = 0.5 - cl.06 
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(f) D = 2.0; Q = 0.5 - c2.06 
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(g) D = 2.6; Q = 0.3 - a3.04 



Time (Myrs) 

(h) D = 2.6; Q = 0.4 - b3.02 



Time (Myrs) 

(i) D = 2.6; Q = 0.5 - c3.09 



Time (Myrs) 



Time (Myrs) 



Time (Myrs) 



(j) D = 3.0; Q = 0.3 - a4.02 



(k) D = 3.0; Q = 0.4 - b4.03 



(1) D = 3.0; Q = 0.5 - c4.03 



Figure 2. Generic examples from the investigated D and Q values. Plots show the evolution of mass segregation (A) with time for = 10 
(solid, includes error bars); = 20 (dashed); N = 50 (dot-dash). The red line indicates a A of unity i.e. no mass segregation. The error 
bars show a la deviation. A movie showing the evolution of (a) al.21 can be downloaded at ^http : / / www. astro . group . shef. ac . uk/stars . ht inT| 
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There are several D = 1.6, Q = 0.3 clusters which ap- 
pear to have little or no mass segregation, for example run 
al.20 (see Fig. [3|. In this cluster there appears to be no sig- 
nificant mass segregation at any time during the 4 Myr of 
the simulation for A'^ — 10, 20 and 50. Fig. 3(b) shows the 
evolution of A for run al.20 in more detail than Fig. 3(a) In 
this plot we show A for A'' = 4, 6, 8, 10 and 12. It is now clear 
that the cluster does mass segregate, but only for stars more 
massive than the G**" most massive star. In fact, a substan- 
tial multiple system consisting of the 4 most massive stars 
is formed in this cluster. Mass segregation is only present 
for the 6*'' most massive stars because the T"' most massive 
star is ejected early in the simulation, thereby enlarging the 
length of all MSTs that include more than 7 stars. Mass 
segregation for = 4 occurs at 0.8 Myr with A — 3.4tj g, 
and by 2 Myr A for the four most massive stars has risen to 
SO.glig.g, with a total MST length of 13000 AU (an average 
separation of 3250 AU), A is very high here because many 
lower mass stars have been widely spread from the cluster, 
making the random MSTs large. This also explains the large 



error values. Fig. 3(b) illustrates the rapid variations in A 



for the small- A*' system of the most massive stars, and the 
fairly smooth variation for larger A'. 

Superficially, run c4.10 looks very similar to run al.20, 
for A'^ — 10, 20 and 50, in that it shows no evidence for mass 
segregation throughout the simulation. When we investigate 
this cluster in detail, we find that the two simulations are 
completely different as this cluster does not undergo a dense 
collapse phase. It becomes only slightly mass segregated at 
3.2 Myr with A for A^ = 4 of S.ltH- 

In runs al.08 and al.21 the stochasticity of runs with 
different random number seeds is easily seen. Both simula- 
tions are initially statistically the same, but in Fig.|3]we show 
the states of the clusters at ~ 1.9 Myr. Run al.08 shows no 
evidence of mass segregation, whilst run al.21 shows very 
significant mass segregation out to A'' = 20 — 25 (in Fig. 4(b) 



the 10 most massive stars are very closely clustered in the 
centre, and the symbols indicating their positions overlap 
somewhat). 

It is important to note that the plots used in Fig. [2] do 
not show all of the mass segregation information. They only 
show mass segregation for the A'^ which is plotted (in this 
case A' = 10, 20 and 50). This choice is rather arbitrary and 
fails to provide all of the information. We have chosen to 
plot these values for clarity, but any detailed investigation 
of the evolution of cool, clumpy clusters requires an analysis 
of a three dimensional dataset including time. A'', and A. 

In run al.04, the stochastic and transient nature of mass 



segregation can be seen. Fig. 5(a) shows that at 2.1 Myr the 
twelve most massive stars in the cluster are in a heavily 



mass segregated state. However, Fig. 5(b) shows that only 
0.1 Myr later the amount of observed mass segregation in 
the cluster has been reduced from A ~ 21 to A ~ 6 because 
of the disruption of the massive multiple system at the core. 

As observations are a snapshot of one time in a clus- 
ter's evolution, they only provide two dimensions of A'^ and 
A. This means that observations miss the time evolution of 
the system which, as is clear in many plots in the supplemen- 
tary material, can vary extremely rapidly. Hence we stress 
that observations only provide a snapshot of a young clus- 
ter. Young clusters can evolve very rapidly on a timescale 
of ~ 1 Myr. Therefore the instantaneous state of a clus- 



1 2 3 

Time (Myrs) 



(a) 




Time (Myrs) 



(b) 



Figure 3. Evolution of A with time for run al.20. Top: A'^ = 
10, 20 and 50. Bottom: N = A (solid line), 6, 8, 10, and 12 (top to 
bottom) for run al.20. Error bars have only been included for 
N = A for clarity. 



ter (especially if it is based on two dimensional positions) 
may provide no information on what the cluster was like 
in the recent past, or to how it mi ght evol ve in the near fu- 
ture (see also 



Iau ison et al.l 



Goodwin fc BastianlboOR : .Bastian et al.ll2008l : 
2009bh . 



3.4 Post-collapse evolution 

The collapse of clumpy clusters may not lead only to mass 
segregation, but also to other interesting phenomena; most 
notably the formation of 'Trapezium-like' systems, the rapid 
ejection of massive stars, and the disruption of the cluster 
itself. 



3.4-1 Formation of 'Trapezium-like' systems 

The Trapezium in the Orion Nebula Cluster (ONC) is a mul- 
tiple system consisting of four of the most massive stars in 
the ONC. For the following discussion we shall use a sim- 
ple definition of a 'Trapezium-like' system as a close group 
of massive stars (at least 3 with M > 8 M©) which are 
bound. Under this criterion, 45 of the 50 -D = 1.6, Q = 0.3 
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Figure 4. Runs al.08 (a) and al.21 (b): D = 1.6, Q = 0.3. The 
ten most massive stars are depicted by triangles. The most mas- 
sive stars are not always clearly seen because of close grouping. 
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Figure 5. Run al.04: D = 1.6, Q = 0.3. The ten most massive 
stars are depicted by triangles, (a) At 2.1 Myr the cluster shows 
very high levels of mass segregation for the twelve most massive 
stars (A 21). The ten most massive stars are tightly grouped in 
the centre of the the cluster, and cannot be seen clearly, (b) 0.1 
Myr later, the amount of mass segregation present in the cluster 
has been vastly reduced (A ?s 6). There are two massive star- 
massive star binaries in the central region, causing symbols to 
overlap. 



Figure 6. Run al.02: D = 1.6, Q = 0.3. The ten most massive 
stars are depicted by triangles, (a) The cluster has evolved to 
form a centrally concentrated distribution. In the core there is 
a tight triple system consisting of 3 of the most massive stars, 
which cannot clearly be seen, (b) 0.7 Myr after (a), the cluster 
has been destroyed by the decay of the multiple system, and most 
of the 10 most massive stars have been ejected. 



simulations formed a 'Trapezium-like' system. These simu- 
lations contained no primordial binaries, but were able to 
migrate the most massive stars into very close proximity 
to each other on very sho rt timescales (~ 1 M yr) through 
dynamics alone (see also, I Allison et al.l [2009bl ). In the 50 
D = 1.6, Q = 0.3 simulations, the average shortest sepa- 
ration reached by the four most massive stars was 0.04 pc 
(about 8000 AU), which occurs within the first 2 Myr. The 
shortest separation between the four most massive stars in 
any simulation was 0.0025 pc (500 AU). We would expect 
that with different initial conditions (i.e., initial binaries, 
smaller Ro, etc) these multiple systems could be made even 
more extreme. In fact, preliminary investigations of simula- 
tions with an initial fractal radius of 0.5 pc (half the size used 
in this paper) and only single stars found a 'Trapezium-like' 
system with the four most massive stars in a cluster with an 
average separation of only 20 AU between them. 

'Trapezium-like' systems form during the dense core 
phase as the most massive stars loose kinetic energy through 
two-body encounters with lower-mass stars. As they mass 
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segregate, the typical separations between tiie most mas- 
sive stars decrease and they form high-order multiple sys- 
tems. We note that the dense core that forms at the peak 
of the collapse is similar to the ini tial conditions used by 
IPflamm-Altenburg fc Kroupal l|2006l ') to model the formation 
of the Trapezium in the ONC. It also has the e xtreme density 
requi red to explain o bservations of ri Ch am (|Moraux et al.l 
120071 1 and the ONC (|Parker et al.ll2009^ . We will examine 
the formation and evolution of 'Trapezium-like' systems in 
more detail in a subsequent paper (R. J. Allison et al., in 
prep.). Hydrodynamical simulations also show the forma- 
tion of trapezium-like sy stems in the gas-dominated phase 
of clu ster formation (see, iBonnell et allbOOSl : iKlessen et al.l 
I2OO9I ). These systems are able to form due to the ability 
of gas to dissipate kinetic energy and redistribute angular 
momentum. Therefore, if the effects of gas were included in 
these simulations we would expect the trapezium systems 
which form in these simulations to appear earlier and to be 
even more extreme. 



3.4.2 Core collapse? 

Some simulations appear to show core collapse events. Early 
mass segregation can create a core-halo structure in the clus- 
ter, with a dense core of massive stars, surrounded by a halo 
of lower-mass stars. The dense core often forms a Trapezium- 
like system (see above). The lower-mass stars are typically 
on radial orbits, and may enter and extract energy from the 
core causing it to shrink and the higher-order massive multi- 
ples to harden. This also causes the halo to heat up and the 
cluster to start dissolving. At some point, the massive mul- 
tiples may decay, ejecting high-mass stars at high velocity 
(see below) and blowing the cluster apart. 

In run al.02, for example, the cluster collapses and 
forms a dense core. The 4 most massive stars form a 
Trapezium-like system which very rapidly decays and de- 
stroys the cluster. Fig. 6(a) shows the cluster during the 



'core phase', the massive star multiple system is in the cen- 
tre of the cluster (but the individual stars are too close to 
resolve in this image). Fig. 6(b) shows the cluster 0.7 Myr 



later. Most of the massive stars have been ejected from the 
central 1 pc, and the cluster has been destroyed. Prior to 
dissolution, the degree of mass segregation increases as the 
central multiple hardens, and at around 1.9 Myr the central 
multiple finally decays, ejecting the most massive stars and 
creating a hard binary. This system contains two stars of 
mass 45.3 and 27.8 M0, which had an initial separation of 
800 AU before the ejection, and a separation of ~ 100 AU 
after. This binary increases its binding energy by around 
1 X 10^° J, which is comparable to the entire potential en- 
ergy of the cluster of ~ —6 x 10^® J. In this simulation the 
formation of the binary system alone is able to directly dis- 
rupt the cluster. It is also possible that the energy input 
comes from gradual interactions with massive star multiples 
which slowly disrupt the cluster and form a hard massive 
star binary. 



Table [5] shows the number of simulations for which the 
final (4 Myr) virial ratio is greater than unity (i.e., the clus- 
ter is unbound). Energy is conserved in the simulations, but 
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Table 2. Fraction of clusters which become globally unbound 
with varying initial virial ratio, Q, and fractal dimension, D. 
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Figure 7. Distribution of final speeds of all stars from the 50 
D = 1.6, Q = 0.3 simulations. The black region shows stars with 
M > 10 Mq , and the dashed line shows the average escape speed 
('^esc ~ 2.5 km s~^). 



an unbound cluster may be formed from an initially sub- 
virial cluster through the redistribution of energy in the core 
collapse phase. 

Even if the cluster manages to survive for the 4 Myr of 
our simulations, shortly after this the most massive stars will 
become supernovae and the loss of the most bound portions 
of the cluster should destroy the cluster. Thus clusters can 
be destroyed dynamically without the need for gas expulsion 
(|Goodwin fc Bastian|[200^ : lGoodwin|[2009l ). 



3.4.3 Massive star ejections 

The ejection of massive stars from clusters can occur in three 
ways. Firstly, Trapezium-like syst ems are inherently unsta- 
ble and can decay dynamically jSterzik fc DurisenI Il998l ) . 
Secondly, other massive stars can interact with Trapezium- 
like systems or with massive binaries and be ejected. Thirdly, 
the hardening of Trapezium-like systems in core collapse 
processes may decrease its stability and cause it to decay 
faster than it might otherwise have done. 

Figure [7] shows the final (at 4 Myr) velocity distribution 
of the 50 D — 1.6, Q — 0.3 simulations, for comparison the 
typical velocity dispersion of a cluster is ~ 1.9 km s~^, and 
the typical escape velocity is ~ 2.5 km s~^. The plot shows 
that it is quite possible for the massive star ejections to reach 
velocities in excess of 10 km s~^. The simulations produce 
stars with masses > 30 Mq having ejection velocities > 10 
km s~^, and in one case a 12 MQStar is ejected with a veloc- 
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ity of 25.9 km s ^ . In all, 13 per cent of stars > 10 M© were 
ejected with velocities > 5 km s~^(~ 5 pc Myr~^). 



4 DISCUSSION 

Many young clusters are observed to be 'mass segregated', 
that is their most massive stars are concentrated pref- 
erentially towards their centres. The ONC (2-3 Myr) is 
well k nown to be mass segregated. iHillenbrand fc HartmannI 
find, using the cumulative distribution of different 
mass groups, evidence of mass segregation for stars more 
massive than 5 M© , and a possibili ty that mass segrega- 
tion is present for stars >l-2 Mq. lAllison et al.l (|2009al ') 
show, using the minimum spanning tree method, that mass 
segregation in the ONC is split into three 'levels'. The 
first contains the four most massive stars, the second stars 
> 5 M0 , and the third is a decrease towards no mass segre- 
gation below 5 Mq. U sing the radial variation o f the IMF, 
iHaravama et all |2008l ') and lSharma et all ()2008l ) show that 
NGC 3603 (2.5 Myr) and NGC 1893 (4 Myr) are mass seg- 
regated, but that there is no eviden ce for mass segregation 
for sta rs < 1 — 2 M© in either cluster. iRaboud fc MermilliodI 
l|l998h show, using the cumulative distribution method, that 
NGC 6231 (4 Myr) is mass segregated. This cluster, like the 
ONC, has different 'levels' of mass segregation. The most 
massive stars (> 16 M©) are clearly more mass segregated 
than stars of lower mass, and stars with masses between 
2.5 and 16 M© are spatially well mixed, but more mass seg- 
regated than stars < 2.5 Mq . Also u sing the cumulative dis- 
tribution method IJose et al] (|2008h find evidence for mass 
segregation in Stock 8 (1-5 Myr), and show th at it is only 
appar ent for stars > 1 M0 . Using our method, ISana et al.l 
l|2010ll find that the cluster Trumpler 14 (~4 Myr) also 
shows mass segregation. There is significant mass segrega- 
tion for stars more massive than «10 M©, but no evidence 
for stars less massive than this. 

Such observations are a strong indication that these 
clusters have undergone an early dense phase. The mass 
segregation in these clusters follows a similar profile - mass 
segregation is present for the most massive stars, and below 
some particular mass ceases to be observed - this is what 
we observe in our simulations. 

There are two important elements that are missing from 
our simulations: binaries and gas. Binaries have been ne- 
glected as we wish to explore the basic gravitational physics 
which binaries would confuse. In future papers, we will in- 
troduce primordial binaries into our simulations. Gas has 
been ignored, as it is computationally extremely expensive 
to include and because introducing gas introduces a whole 
new set of parameters (the equation of state of the gas, the 
Mach number and power spectrum of the turbulence, etc.). 
But clearly gas is a vital ingredient in very young star clus- 
ters. It is what the stars form from, and it can be a significant 
contributor to the background potential. 

Prom our simulations, it is not clear if mas- 
sive stars might fo rm directly in massive cores (e.g., 
iKrumholz et al .ll2007t). through fragment-induced starvation 
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2010) or due to competitive accretion (e.g., 
ISmith et al.l 2009 1. However, hydrodynamical simulations 
indicate that gravoturbulent fragmen tation leads to sub- 
clustered and mass segregated clusters (|Mac Low &: KlessenI 



I2OO4I : iMaschberger et al.ll20ld ) . Clearly, the ability of young 
clusters to dynamically mass segregate on a very short 
timescale means that massive stars do not have to form (or 
rather, build up their mass) in the core of a young cluster. 
Massive stars can form in subclusters or in the outskirts of 
clusters and within a Myr or so be very mass segregated. 
However, not only stars, but also gas, will be involved in the 
collapse of the cluster to a dense phase and so cool collapse 
might enhance competitive accretion. Indeed, it is difficult 
to imagine that some competitive accretion will not occur 
during the channelling of stars and gas into a very dense 
state. Do massive stars form in massive cores that are then 
dynamically mass segregated? Does competitive accretion 
dominate, especially during the dense collapse to form the 
massive stars? Do fairly massive stars form in the outskirts 
of clusters which are dynamically mass segregated and then 
have their masses increased by competitive accretion in a 
hybrid scenario (such a process may be occurring in simula- 
tions; see Maschberger et al. 2010)? 



5 CONCLUSIONS 

Observations and theory both strongly suggest that the ini- 
tial conditions of star clusters are cool and clumpy. There- 
fore, we have conducted a large number of A''-body simula- 
tions of the early (< 4 Myr) evolution of clusters with virial 
ratios of Q = 0.3,0.4 and 0.5 (where 0.5 is virialised), and 
fractal dimensions of D = 1.6,2.0,2.6 and 3.0 (where 3.0 
is roughly a uniform density sphere) with radii of 1 pc and 
1000 members (ie. total masses of ~ 500 M©). In our simu- 
lations, all members were initially single stars selected from 
a Kroupa IMF, and the simulations lasted for 4 Myr (so not 
requiring stellar evolution to be included). 

We study the evolution of the star clusters with a partic- 
ular emphasis on the level of mass segregation. We measure 
mass segregation using a minimum spanning tree to provide 
a quantitative measure of mass segregation and which is not 
biased by clumpy underlying mass distributions (Allison et 
al. 2009a). 

This study follows that of Allison et al. (2009b), in 
which we showed that clusters with Q = 0.3 and D = 1.6 
(i.e., very cool and extremely clumpy) undergo collapse to 
a short-lived but extremely dense core. This core can dy- 
namically mass segregate the most massive stars down to a 
few Mq via two-body encounters before re-expanding. The 
re-expansion is partially driven by the increase in the ve- 
locity dispersion of the low-mass stars caused by two-body 
encounters. 

Our main results may be summarised as follows: 

• The depth of the collapse, and so the degree of mass 
segregation depend on both the initial virial ratio, Q, and 
the degree of substructure, D. Low-Q and low-D clusters 
can collapse to a denser state and so mass segregate more 
than high-Q, high-D clusters. 

• Whilst there is a general trend of increasing mass seg- 
regation with lower-Q and lower-D the inherently stochastic 
nature of fractals means that statistically identical clusters 
may undergo very different evolution. 

• Many features are extremely short-lived, and young 
clusters can change rapidly and violently. Observations only 
provide a snapshot of the evolution of a cluster and it is 
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dangerous to draw conclusions about the past and future 
state of a cluster from a single snapshot (see also Bastian et 
al. 2008). 

• Young clusters can undergo core collapse. Early dy- 
namical mass segregation establishes a core of massive 
stars (see also Pflamm-Altenburg & Kroupa 2006) and a 
halo of lower-mass stars. Energy loss from the core to the 
halo can drive the formation of massive, hard binaries and 
Trapezium-like multiple systems. The heating of the halo 
and the hardening and subsequent decay of the central mul- 
tiple systems can dynamically destroy a cluster within a 
few Myr. 

• The interactions of massive stars in the core of a young 
cluster can cause the ejection of even very massive stars at 
velocities in excess of 20 km s~^This may help explain ejec- 
tions from the ONC (see also Pflamm-Altenburg & Kroupa 
2006). 

• The early evolution of cool, clumpy clusters is rapid, 
violent, and extreme. The densities of clusters (and hence 
their crossing and relaxation times) can change by orders of 
magnitude during the first few Myr of their existence. Thus 
the currently observed properties of young clusters are just a 
snapshot in the life of these clusters and extreme care must 
be taken in inferring the past history or future evolution of 
clusters from a single snapshot (see also Goodwin & Bastian 
2006; Bastian et al. 2008; Allison et al. 2009b). 

That young clusters are mass segregated down to a few 
M0 , but not below, is due to the short-lived dynamical mass 
segregation phase which is able to mass segregate only the 
most massive stars. That the ONC has an unstable high- 
order multiple containing four of the most massive stars can 
be explained by its dynamical formation during the dense 
phase. The ejection of high-mass stars from clusters can oc- 
cur during the dense phase, or afterwards from the decay of 
higher-order massive multiple systems. 
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